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ABSTRACT 

We use a new stacking technique to obtain mean mid IR and far IR to far UV 
flux ratios over the rest near-UV/near-IR color-magnitude diagram. We employ 
COMBO-17 redshifts and COMBO-17 optical, GALEX far and near UV, Spitzer 
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IRAC and MIPS Mid IR photometry. This technique permits us to probe infrared 
excess (IRX), the ratio of far IR to far UV luminosity, and specific star formation 
rate (SSFR) and their co-evolution over two orders of magnitude of stellar mass 
and redshift 0.1 < z < 1.2. We find that the SSFR and the characteristic mass 
(Mo) above which the SSFR drops increase with redshift (downsizing). At any 
given epoch, IRX is an increasing function of mass up to Mq. Above this mass 
IRX falls, suggesting gas exhaustion. In a given mass bin below Mq IRX increases 
with time in a fashion consistent with enrichment. We interpret these trends 
using a simple model with a Schmidt-Kennicutt law and extinction that tracks 
gas density and enrichment. We find that the average IRX and SSFR follows 
a galaxy age parameter ^ which is determined mainly by the galaxy mass and 
time since formation. We conclude that blue sequence galaxies have properties 
which show simple, systematic trends with mass and time such as the steady 
build-up of heavy elements in the interstellar media of evolving galaxies and the 
exhaustion of gas in galaxies that are evolving off the blue sequence. The IRX 
represents a tool for selecting galaxies at various stages of evolution. 

Subject headings: galaxies: evolution — ultraviolet: galaxies 



1. Introduction 



It has long been recognized that the present day pro perties of mo s t galaxies can be 
repre sented by relatively simple star formation histories (ITinsleyl Il968l : ISearle fc Sargent 
19721 ). While the physical basis for exponential star formation histories is almost certainly 
oversimplified, the resulting spectral energy distributions generally do an excellent job in 
representing galaxy spectra and broadband colors. If exponential models have a basis in the 
physics of star formation history, in particular in the conversion of gas into stars, then they 
make basic predictions that relate the specific star formation rate to the gas fr ac tion in galax- 



i es oy er time. Evidence for such a picture is growing (e.g., iBell et al.l (120051 ): iNoeske et al. 
( 120071 )). Coupled with a model for chemical evolution, and brushing aside for the moment 
the complexities of dust reprocessing, this framework could also provide a description of the 
evolution of dust extinction in galaxies and the co-evolution of extinction and star formation 
rate. In particular, we could expect a growth in the dust-to-gas ratio as gas is processed 
through star s and potentially an increase in extinction over time (e.g., as seen at z~2 by 
Reddy et al.l (120061 )). At the same time, as galaxies exhaust their gas supply, by whatever 



mechanism, we may detect a corresponding drop in extinction. 



In order to discern such effects, we need to segregate galaxies by a parameter which is 
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likely to be related to the timescale for evolution. There is certainly theoretical motivation 
for using stel lar mass as the fundame ntal parameter. For example, surface density scales with 
stellar mass jKauffmann et al. 2003 ). and star formation rate scales with gas surface density 
(lKennicuttiri989h. T he ob servational case for "downsizing" s eems secure (ICowie et al.lll996l 



Brinchmann fc Ellis! 120001 ). The mass-metallicity relation (ITremonti et al.l l2004j ) suggests 
that low metallicity in lower mass galaxies could be related to higher gas fractions and 
lower processing through sta r formation. Lower mass galaxies have younger stellar ages 



Kauffmann et al.l (120031 )). Ther e is a well-known relation ship between luminosity as a 



(e.g. 

proxy for stellar mass and extinction (IWang fc Heckmanlll996f) which is present even at high 



2000 



Reddy et al.ll2006l : iPapovich et al. 



redsh ift ( iMeurer et al.l Il999l : lAdelberger &: Steidel 
20061 ). We have recently established a tight relationship between m etallicity and infrare d 



excess (IRX), the ratio of far infrared to far ultraviolet luminosity (j Johnson et al.l 120071 ). 
that suggests that IRX may be used as a tracer of metallicity and its evolution. Finally, 
there is growing evidence that the so-called "blue cloud" of star form ing galaxies on the 



color-magnitude diag ram is actually a "blue sequence" in stellar mass (IWyder et al.l 12006 



Johnson et al.ll2006al ) that is relatively tight in color space when extinction is corrected. 



In order to distinguish trends with stellar mass, it is critical to have as large a mass 
dynamic range as possible. At the same time, dust extinction is likely to be a complex 
process that introduces considerable noise into any overall trends. Inclination variations 
alone produce much variance for an otherwise constant dust geometry and extinction law. 
We need to develop an approach which reveals the average trends with stellar mass in spite 
of this noise. A major benefit of a large multiwavelength survey is the ability to extract 
such trends by averaging over many galaxies. We have used Spitzer IRAC data to measure 
stellar mass and MIPS24 data for dust luminosity. We combine this with GALEX UV and 
COMBO-17 optical photometry and redshifts. A major difficulty we face when combining 
IR and UV survey data is the relatively small overlap in detected sources, with the bulk of 
the overlap occuring at high luminosity and mass. Thus we have developed a new stacking 
approach which permits us to study IRX over two orders of magnitude in stellar mass over 
the redshift range 0.05 < z < 1.2. Using this and a bolometric correction we obtain an 
average IRX over the UV, H-band color magnitude diagram. We use this to find total star 
formation rate, stellar mass, and specific star formation rate. Finally, we show that the co- 
evolution of the average IRX and specific star formation rate can be modelled using simple 
exponential star formation histories and closed box chemical evolution to ^ ~ 1. 



We note that 



Zheng et al.l (120071 ) have recently used an independent stacking technique 



( IZheng et al.l 120061 ) to derive the SF history vs. stellar mass, also using COMBO-17 and 
Spitzer data, and reached many conclusions that are similar to ours, although with important 
differences which we discuss in 55.2. 
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We use a concordance cosmology = 0.70, fi^ = 0.30, and Hq = 70 km s~^ Mpc^^ 
We use AB magnitudes for all bands. We also use the following nomenclature: observed 
magnitudes are given by rrii, for example the observed NUV and R magnitudes are rriNuv 
and rrifi. Rest-frame magnitudes are denoted by FUV, NUV, H, etc. Extinction corrected 
rest-frame absolute magnitudes are denoted Mfuv,o, ^nuv,o, ^h,o, etc. Finally, we define 
the infrared excess IRX as the log ratio of FIR to FUV luminosity iyCy), unless specifically 
called out. 



2. Data and Source Catalogs 



2.1. Primary Datasets 

2.1.1. GALEX 

The GALEX observations of CDFS consist of a total of 61 orbital visits over the pe- 
riod 4 November 2003 to 5 November 2005 for a total exposure time of 49,758 seconds. 
Simultaneous exposures were obtained in the Far UV (FUV, 1344-1786A, center 1549A) 
and Near UV (NUV, 1771-283lA, center 2316A) bands. The individual and co-added im- 
ages were processed using version 5.0 of the GALEX data pipeline, also used to process 
GALEX data releases GR2 and GR3. The 1.25 ° diameter GALEX images completely cir- 
cumscribe the other two survey footprints. The GALEX mission, on-or bit performance, an d 
current status o f calibr ation and pipeline reductions are summarized in iMartin et al.l (120051 ) , 
Morrissey et al.l (120051 ) . and iMorrissey et al.l (120071 ) respectively. Source photometry errors 
(systematic) should be less than 0.05 mag and astrometric errors less than 1 ". Images of 
this exposure level in low background regions should reach a 5cr Poisson limited depth of 
rriNUV —25.5 AB magnitudes in both bands, and roughly 3a at nriNuv — 26.0. However NUV 
data in particular are confusion limited because of the 5-6 " PSF FWHM. We therefore used 
a PSF fitting source extraction procedure that uses the CDFS optical positions as priors. 
This is described below. 



2.1.2. COMBO- 11 



The COMBO-17 survey (jWolf et al.l l2003l ) combines a set of medium and wide photo- 
metric bands to obtain robust photometric redshifts and basic object classification to a depth 
of ~24. A complete description of the survey can be found in I Wolf et al.l (120041 ). The 
CDFS field is 0.5 by 0.5 degree^ centered on (a,5)j2ooo = (03''32'"25^ -27°48'50"). Other 
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than redshifts, we use the COMBO-17 survey for two purposes: in order to generate a k- 
corrected NUV luminosity and NU V-H color, and as the basis for the PSF fitting extraction 
of the FUV and NUV source fluxes. IWolf et al.l (120041 ) have used a Monte-Carlo technique to 
derive the survey completeness vs. object color, type, and magnitude. We have used these 
completeness matrices to derive the volume-corrected distributions, as we describe below. 



2.1.3. Spitzer 



The Spitzer data are described in detail in iPerez- Gonzalez et al.l (120051 ). which we 
briefly recap here. The 1°.5 x 0°.5 rectangular areas centered on CDFS, (a,5)j2ooo = 
(03''32'"02^ -27°37'24"), are mapped in MIPS 24 //m in scan map mode, and also in the four 
IRAC channels (3.6, 4.5, 5.8, a nd 8.0 iim). The MI PS 24 fim reduction was performed using 
the MIPS Data Analysis Tool (IGordon et al.ll2005l ). resulting in images with average expo- 
sures of ~1400 s. IRAC images were reduced with the general Spitzer pipeline and mosaiced, 
yielding an average exposure time of 500 s. Source catalogs for IRAC 3.6 /im detections are 
used below i n a jointly selected san iple. We tested catalogs generated by a simple one-pass 
SExtractor (IBertinfc Arnoutd Il996l ) extraction, and by the multiband technique used by 
Perez-Gonzalez et al. J2005h . with no significant differences noted in our results. Source cat- 
alogs of MIPS 24 fim objects were used to clean 24 /xm images for stacking, as we describe 
below in §3.3[ Again, a single-pass SExtrac tor catalog produced very si milar results to the 
mulit-pass PSF-fitting catalog generated by iPerez-Gonzalez et al.l (120051 ). 



2.2. Matched Datasets 

2.2.1. GALEX/COMBO-17 PSF Fitting Catalog 

As noted above, deep GALEX images suffer from source confusion, especially in the 
NUV. For fields with complementary deep optical photometry, we can use the positions of 
sources from the optical catalog to deblend the GALEX images and obtain more reliable 
flux estimates. The center of the COMBO-17 field is only slightly offset (3.8 arcmin) from 
the center of the GALEX images and is much smaller than the GALEX image. Within the 
COMBO-17 field, the variation of the GALEX PSF is small, and so we have have used one 
average PSF for each band. After correcting for the small (less than 1 arcsec) systematic 
offsets between the GALEX and COMBO-17 astrometry, the deblending proceeds by dividing 
the region to be deblended into contiguous 100 x 100 pixel chunks and then simultaneously 
fitting the amplitudes of the sources at positions taken from the optical catalog and the 
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mean background in each chunk. We assume that the counts in each pixel are Gaussian- 
distributed, which is a safe assumption for the NUV (where the background level is 100 
counts), but is questionable for the FUV (where the background level is 10 counts). In 
order to test the reliability of our deblending, we have added approximately 1000 artificial 
point sources to the C0MB0-17/GALEX overlap region and then compared the extracted 
fluxes with known input fluxes. In the FUV, the deblended magnitudes are systematically 
fainter than the input magnitudes by 0.04 mag and have errors that are 20% larger than 
expected from counting statistics. In the NUV, the deblended magnitudes are systematically 
too faint by only 0.01 mag, but the errors, due to the source crowding, are a factor of 2 larger 
than expected from counting statistics. For the 49,758 second GALEX images used here, 
the 95% confidence detection limits are 25.55 mag in FUV and 25.10 mag in NUV. 



2.2.2. Merged Catalog 

There is a low fraction of sources detected in all three catalogs. We generated individual 
source catalogs for the five Spitzer images using Sextractor. We matched these detections 
to the merged C0MB0-17/GALEX catalog using a 2"search radius. The common area of 
the three surveys is 0.19 degree^. The matched source statistics are summarized in Table 
[U Of the 11778 COMBO-17 tjir <24 sources in the common region, 7498 are detected at 
mNuv <26, 1784 in MIPS24 (>0.02 mJy), 4955 in NUV and IRACl (>0.0002 mJy), and 
1171 in NUV, IRACl and MIPS24. As we discuss below, the main explanation for the 
low overlap fraction is that many UV-selected sources have moderate to low infrared-to-UV 
ratios and are not directly detected in the mid-infrared. 

Because we are keenly interested in the evolution of the average extinction, infrared 
excess, and star formation history of galaxies over cosmic time, we have adopted a stacking 
approach. We summarize the complete methodology in the next section. 



3. Analysis 



Our goal in this paper is to determine the evolution of the average IRX and extinction 
and relate this to the evolution of the star formation rate, as a function of stellar mass. We 
would like to exploit a property of the blue sequence of star forming galaxies that is rapidly 
becoming clea r, that this sequen c e has a relatively l ow dispersion of pr operties once the 
mass is given (INoeske et al.l 120071 : IWyder et al.l l2006l : iMartin et al.l l2007al ) . The dispersion 
of various properties, such as extinction, stellar age and mass, measured in individual bins of 
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the NUV-r color magnitude diagram, for example, is low. The IRACl channel and COMBO- 
17 R-band give a good estimate of the rest H-band flux over the redshift range < z < 1.2, 
providing a good stellar mass tracer with low extinction sensitivity. We therefore feel it is 
reasonable to stack using bins of the rest-frame M^, NUV-H band color- magnitude diagram. 



3.1. Methodology 

Here is a step-by-step summary of our approach, with cross-references to more detailed 
discussion: 

1. Use COMBO-17 positions to generate a joint C0MB0-17/GALEX catalog using PSF 
fitting (gMH) 

2. Match IRACl sources detected with Sextractor to joint C0MB0-17/GALEX catalog 

(gIZ2D 

3. Generate rest-frame NUV, H-band magnitudes using SED interpolation ( §3.2p 

4. Construct the volume-corrected color-magnitude diagram (CMD) 0(Mh ,NUV-H) in 
several redshift bins 



5. For each (Mr , NUV-H, redshift) bin, stack all IRAC2-4 and MIPS24 images at the 
R-band COMBO-17 source positions falling in that (Mh , NUV-H, z) bin. Our stacking 
technique adds detected source fluxes to a stack of undetected source regions, as we 
discuss in ^3.31 



6. Use a bolometric cor rection obtained by fitting a local SWIRE/SDSS/GALEX sample 



(jjohnson et al.ll2007l ) and the total Far IR luminosity and infrared excess lRX=log Lfir/ Lpuv 



or log Lfir/Lnuv- {^M 

7. Use a standard extinction law to convert IRX into NUV and H-band extinction. ( §3.41) 

8. Determine the average extinction correction for galaxies in each (Mh , NUV-H, z) bin. 

9. Using {NUV — H)o and Mh^ infer the specific star formation rate and stellar mass 
using a simple prescription. ( §3.41) 

10. Determine the volume-corrected distributions of specific star formation and stellar mass 
(M) in each redshift bin 0(A^,SSFR,z). f TOj) 
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11. Finally, calculate the average IRX and specific star formation rate {{IRX) and 
(SSFR)) as a function of stellar mass and redshift ( §4.31 14.41) . 



3.2. Evolution of the Color-Magnitude Distribution 

We choose to use rest-frame FUV or NUV to derive star formation rate and H-band 
magnitudes to obtain stellar mass. Rest-frame H-band flux was obtained by interpolating 
between the COMBO-17 R-magnitude and the IRAC 3.6 fim flux, exploiting the fact that the 
SED is essentially constant over this range for most galaxy templates. Rest-frame NUV flux 
is an interpolation of the observed NUV and the COMBO-17 catalog rest-frame u-band flux. 
Rest-frame FUV flux is an interpolation between observed FUV and NUV which accounts 
for the Lyman continuum break. We also tried SED fitting, which produces very similar 
results. 

We derived the volume-corrected, rest-frame NUV-H vs Mh distributions as follows. We 
used five redshift bins, z =0.05-0.2. 0.2-0.4, 0.4-0.6, 0.6-0.8, and 0.8-1.2. Maximum detection 
volumes in were derived for a sample jointly selected in observed NUV, r-band, and IRAC 
channel 1 with the following limits: niNuv < 26.0, rrir < 24.0, and /(3.6) > 0.5/iJy. The 
minimum Vmax of the three bands was used. For each galaxy a completeness was calculated 
fc = fc{NUV) * fc{C17). We assume that the IRAC channel 1 completeness is high to the 
flux limit. With PSF fitting, NUV completeness is estimated to be ~80% at m7vc/F=25.5 and 
~56% at mAr[/y=26.0, based on comparison between the observed magnitude distribution 



and that measured by I Gardner et al.l (l2000l ). Because of the soft rolloff of completeness 



for this PSF-fit catalog we choose to use a deeper magnitude cutoff corresponding to a 3a 
detection threshold. The COMBO-17 redshift catalog completeness is a function of object 



magnitude, color, and type. We used the completeness matrix derived by IWolf et al.l (120031 ) 
to calculate /c(C17) for each galaxy. The volume-corrected color-magnitude distribution is 
then calculated by summing for each galaxy the term 



Vmax,i = min[Vmax{NUVi, Zi), Vmaxiru Zi), Vmax{f{3-G)i, Zi)] fc{NUVi) fciVi] Qi - Ti) [l] 



<P{NUV - H, Mh, z) = Y^ tf^ (2) 

''maxA 



The resulting distribution is displayed in contour plots in Figures [H The general trend 
that can be seen is a shift to bluer NUV-H colors and brighter Mh magnitudes. These plots 
also show average IRX in each bin, to be discussed in the next section. 



- 9 - 



3.3. IR Stacking, Bolometric Correction, Extinction Correction 



We have generated an average IRX for each bin in the Mh , NUV-H color magnitude 
diagram. As we discussed earlier, we do this because of the small overlap in UV and MIPS 
detections. There is considerable evidence that galaxies occupying a single col or magnitude 



bin have a relatively smal l dispersion in most pr operties, including extinction (IMartin et al. 



2007al : IWyder et al.ll2006l : I Johnson et al.ll2006al ). This approach allows us to estimate the 
average IRX over a large range of stellar mass and redshift, offering sensitivity to quite low 
IRX. In order to ensure that the technique is not affected by systematic or random error, we 
perform a set of tests below. 

The basic stacking technique is as follows. We first generate a catalog of detected 
sources in the MIPS 24 fim band using Sextractor. Using this catalog we generate a set of 
cleaned images with detected sources removed. For each redshift and each color-magnitude 
bin (NUV-H, Mh , z) we stack images in each band that do not have detected sources. We 
then extract either a source flux or an upper limit, and add this to the detected flux. This 
results in a flux or upper limit. 

We must then make a bolometric correction to the observed 24 /im l umin osity. We have 
used the GALEX/SWIRE catalog generated for Ijohnson et al.l (j2006al . 120071 ) to derive the 
bolometric correction of rest frame flux from 12-24 /xm (corresponding to 0< z <1). We use 
the measured FIR fluxes and flts to iDalefc Heloul (120021 ) SEDs and derive coefficients in the 
following relationship using log-log fits: 



log W = ax + bx\og{L[X]/10'"LQ) 



(3) 



where A is the observed rest frame luminosity {L[X\ = z/Lj,[A]). We list in Table 2 the 
coefficients aA and b^. The rms errors in the fits used to derive the coefficients are small, 
a ~ 0.03 - 0.06. 

Finally, we correct the NUV-H color and H-b and magnitude for in ternal extinction using 
the IRX and the following prescription based on ICalzetti et al.l (120001 ) : 



A = 2.51og,o[-^/i?X, + l] (4) 

where i corresponds either to FUV or NUV, and the bolometric corrections are BCdust = 1-75, 
BCfuv = 1-68, and BCnuv = 2.45. 

Since the H-band correction is small, the extinction correction is very insensitive to 
the extinction law. It enters somewhat if we use L(FUV) to generate IKKpuv and use this 
to correct (NUV-H). Even if there is evolution in the extinction law, we have found that 
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using IRXnuv, gives very similar results to those using IRXpuv The volume-corrected 
distribution of extinction-corrected magnitudes 0(Mh 05[NUV-H]o) vs. redshift are given in 
Figure [21 We show in Figure [3] the uncorrected and corrected distribution in NUV-H. 



3.4. Stellar Mass and Specific Star Formation Rate 



We derive stellar mass Ad using extinction correcte d rest-frarne H-b and absolute mag- 
nitude and NUV-H color, following the basic scheme of iBell et al.l (120031 ) . For smooth star 
formation histories the stellar mass-to-light ratio is a single parameter function of a measure 
of the specific star formation rate such as the rest frame extinction corrected NUV-H color. 
This can be seen in Figure HI which shows the predicted A4/C vs. (N UV-H) n for different 



values of the exponential SFR decay, based on solar metallicity models of iBruzual fc Chariot 



(2003) and a (standard, non-diet) Salpeter initial mass functioE0. There is almost no depen- 
dence on the star formation history for (NUV-H)o for NUV-H<2.5, where the bulk of the 
extinction-corrected galaxies fall. We use this parabolic fit: 



\og^o{M/C) = -0.667 + 0.17{NUV - H)o + 0m373{NUV - H)l (5) 

We have also tested more complex star formation histories in which starbursts become 
significant. These models produce the same general trends between M./ C and NUV-H color, 
with some dispersion. There is no significant impact on the results described below. 

We derive the star forma tion rate from th e extinction-corrected FUV luminosity using 
SFR = 1.4 X 10~^^£;^(1500) (jKennicuttlllQQSl ). We obtain the specific star formation rate 



(SSFR) by dividing by the stellar mass. 

We also note that the specific star formation rate (SSFR) is tightly correlated with 
(NUV-H)o and independent of decay timescale, as can be seen in Figure HI We can also use 
the following linear fit to convert (NUV-H)o to SSFR. 

\og^o{SSFR[yr-^]) = -7.8 - O.Q5{NUV - H)o (6) 
Again, either technique for computing the SSFR produces essentially identical results. 

We derive the stellar mass and SSFR for each galaxy by correcting its rest-frame Mh , 



^No te that all derive d stellar masses and star formation rates can be converted to the "Diet Salpeter" 
IMF of iBell et al.l (j2003l ) by multiplying by 0.7 Specific star formation rates are unaffected. 
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Mi500; and NUV-H for the (bin- averaged) extinction. Using these distributions we derive the 
volume-corrected bivariate A^-SSFR distribution (/)(A^,SSFR,z) in the same fashion as the 
volume-corrected color-magnitude distribution. We also generate an average IRX in each bin 
of the A^-SSFR distribution. We calculate the mean log infrared excess in each bin, using 
the IRX obtained in the previous section. These distributions are displayed in Figure [51 



3.5. Errors 



We use the bootstrap method (jEfronlll979l ) to derive errors to all bivariate distributions 



discussed above as well as the average distributions discussed in the next section. Specifically, 
in each redshift bin we randomly select objects, with replacement, until we have the same 
number of objects found in that redshift bin. We then proceed to determine the color- 
magnitude distribution, the mean 24 micron flux in each color-magnitude bin by stacking 
this new sample, the corrected CMD, the Mass-SSFR distribution, and the mean IRX in 
each Mass-SSFR bin (cf. §4.31) . Of order one hundred trials are used to generate a standard 
deviation in each bin of every distribution calculated. Such errors will not, however, account 
for cosmic variance due to large scale structure, (cf. §5.21) . 



Results 



4.1. Infrared Excess vs. Stellar Mass 



We begin by examining the trends in the NUV-H, Mh CMD. At a fixed Mh redder 
galaxies have a higher IRX. In general, the blue sequence shows a significant tilt in the color- 
magnitude diagram, much of which appears to be produced by this extinction-luminosity 
relationship. Much of the color width of the blue sequence, which l eads some authors t o refer 
to it as the "blue cloud", is also produced by vari ance in extinction (jWvder et al.ll2006l ). some 
of which is simply due to inclination variations (IMartin et al.ll2007al ). Extinction correction 
produces a much tighter distribution in the CMD, as we see in Figure [2] and Figure El The 
trend of increasing IRX with H-band luminosity is even more apparent in the extinction- 
corrected CMD. Consequently, there is a strong increase in IRX with stellar mass as is 
expected from the trend in the CMD. This trend can be clearly detected in Figure [51 This 
trend persists in all redshift bins. 
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4.2. Evolution of the Bivariate CMD 



There is clear evolution in the NUV-H, Mh color-magnitude diagram, in the sense that 
the density of H-band luminous galaxies is inc reasing with redshift. Thi s is consistent with 
the increase in characteristic UV luminosity (jSchiminovich et al.l |2005| ) and B luminosity 
( iBell et al.ll2004l ). As expected from the previous section, this is accompanied by an increase 
in the contribution from higher IRX galaxies. The evolutionary trend is even easier to discern 
in the extinction corrected CMD, Figure [2l 



4.3. Infrared Excess vs. Stellar Mass and Redshift 

In order to further explore IRX-mass relationship and its evolution we derive the average 
IRX in each mass and redshift bin. We have calculated this average using the number density 
{IRX){Ai) and weighted by the star formation rate {IRX)sfr{-M.)- The mass trend in 
redshift bins is shown in Figure [71 The average IRX increases sharply with mass up to a 
critical mass. The slope in the in the IRX-log mass relation is greater than one. The critical 
mass is lower at low redshift, ~ 10.5 at z~0.3, but appears to move to higher mass at higher 
redshift, with ~ 11.5 at z~l. 

The redshift trend in mass bins is shown in Figured In the highest mass bin with good 
redshift coverage (logMcrit = 11-5), IRX increases slowly with time then sharply decreases 
for z<0.5. In the lowest mass bin, IRX appears to increase with time to the lowest redshift 
bin. These trends appear in both the number and SFR-weighted average IRX. For our 
subsequent analysis we use the number-weighted average. 



4.4. Co-evolution of SFR and IRX 

We have seen that the star formation rate density is moving to higher masses at higher 
redshift. This can best be seen in the SFR-weighted bivariate A^-SSFR distributions shown 
in Figure [6l It is interesting to determine the average SSFR vs. mass and redshift as we did 
for IRX. The number-averaged SSFR is given in Figure [91 This shows that at lower mass, 
logM. = 9.5, (SSFR) evolves slowly, while at higher masses the SSFR falls rapidly with 
time. 

The behavior in Figures U]M can be explained by a simple model. We suppose that 
average IRX is determined principally by the gas surface density and by the metallicity. 
This naturally produces a rising then falling IRX as the gas becomes enriched (in a closed 
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box model) and then exhausted. Downsizing imphes that star formation, enrichment, and 
ultimate gas exhaustion move to lower masses with time, consistent with these resuUts. 

We examine this model further in the next section. But first we ask whether the observed 
trends could be an artifact of selection effects or other aspects of our technique. 

4.5. Issues and Caveats 

We have performed a number of tests to ensure that the results presented above are not 
a product of the samples or analysis approach. 

We could use either FUV (1530A) or NUV (2270A) flux to derive star formation rates. 
Since we bin and stack sources in the NUV-H CMD, there could be systematic effects intro- 
duced by either the use of an extinction law to correct FUV given the NUV-derived IRX, or 
there could be effects introduced by the different data samples used to derive FUV and NUV 
rest luminosities. The former come completely from GALEX data, while the later come from 
interpolating GALEX and COMBO-17 data. We find however that there is no significant 
difference in the results using FUV or NUV to derive IRX and SFR. 

We tested stacking the MIPS24 data using detected sources and cleaned images, and 
using only the fluxed images (stacking detected and undetected sources together). This 
produced no statistically significant difference. We also were concerned about the MIPS24 
detection limit and whether the low limit of 0.02 mJy used to detect and clean the images 
would include some spuriously detected sources due to confusion, yielding an artificially high 
24 micron flux in the stacked result. We checked this by increasing the detection limit by 
a factor of two and repeating the entire analysis. Again, these results showed only minor 
quantitative changes. 

A very important question is whether our census of objects is complete. We could 
be missing FIR luminous objects that fall below the UV magnitude limits of the sample. 
Moreover, it is likely that as we move to higher redshift, more of the high IRX and/or low 
SSFR sources are lost due to the UV magnitude limit. This could clearly introduce a spurious 
blueing trend as redshift increases, which is exactly what we detect. 

To test for this effect, we repeated the analysis on the following samples: 1. Baseline: 
NUV<26.0, r<24.0; 2. Case 2: NUV<25.0, r<24.0; 3. Case 3: NUV<27.0, r<24.0; 4. Case 
4: NUV<26.0, r<25.0; 5. Case 5: all r<24.0 objects, whether or not detected in NUV. 
Those objects below the NUV detection limit are given an artificial magnitude NUV=27.0. 
The average IRX vs. redshift for these cases is shown in Figure [101 There are no significant 
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changes to (IRX) or in the observed trends with mass and redshift. 

Another test is to consider the inchnation bias of the sample. A sample at high redshift 
which has not included highly inclined, more reddened galaxies of otherwise similar overall 
dust content will display a higher average minor to major axis ratio than the low redshift 
counterpart. The average axis ratio obtained from the (seeing limited) COMBO-17 data 
shows no significant trends with redshift. This is true in the Mh , NUV-H CMD and as we 
show in Figure im the final A^-SSFR diagram. 

Finally, we used a Monte Carlo model to test whether the evolving IRX-mass rela- 
tionship could be an artifact of the sample selection. The model is semiempirical and we 
briefly summarize it here. The model predicts the bivariate luminosity function in the ex- 
tincted, rest-frame NUV-H, Mh CMD, and the distribution of (IRX) over the same CMD 
{{IRX)[Mi{ , NUV-H]). The number distribution is given by a Schechter function in mass. 
The SSFR is log-normal with a constant mean and variance to a certain critical mass, then 
falls. IRX is log-normal and the mean IRX scales with mass. Evolution is introduced into 
the number density, mean SSFR, SFR cutoff mass, and the IRX-mass relationship. In the 
latter case, the following relationship is introduced: 

IRX = IRXo + am (log M/Mo) + a„,, (log M/Mq) log (1 + 2;) + amx * ^ (7) 

where ^ is a normally-distributed random variable. This assumption allows for IRX depen- 
dence on mass, and evolution of this dependence in a mass-dependent fashion, in other words 
the evolutionary trends we appear to detect in the data. 

We convert mass, SSFR, and IRX into observed SED's using (in reverse) the identical 
transformations that we used for the data. SEDs are redshifted and run through detection 
filters with appropriate completeness functions. We then subject the list of objects and 
observed FUV, NUV, R-band, IRAC, and MIPS24 fluxes to the identical analysis steps as 
the actual sources, producing the various distributions. (We do not simulate the actual 
image formation and detection process). Finally, we compare the Monte Carlo and data 
distributions using a chi-square statistic. For this comparison we combine data and Monte 
Carlo errors (data errors calculated from bootstrap, Monte Carlo errors calculated from 
Poisson statistics). To calculate we use all bins in which either data or Monte Carlo results 
are predicted. We simultaneously fit both (j){Mu, NUV - H) and {IRX){Mn,NUV - H) 
over all five redshift bins. 

Because the model has ~ 10 — 14 parameters, it is difficult to guarantee any given 
local minimum is the global minimum. Extensive experimentation has shown that the 
distribution and IRX distributions are mainly influenced by separate variables, so some 
minimization can be decoupled. We find best fits with 0^,2 —1.5 and — 0.3, with 
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log Mo ~ 11.0 and with small formal errors (< 0.05). The latter are derived in the usual way 
by fixing the parameter of interest and marginalizing over all others. The key conclusion 
is that the significant non-zero value of am,z provides additional evidence that the evolving 
IRX-mass relationship is not an artifact of the sample selection. 



5. Discussion 

5.1. Simple Extinction, Metallicity, and Star Formation Evolution 

The phenomena displayed in Figures [8] and [7] have a very simple interpretation. At 
low mass, ongoing enrichment by star formation increases the dust-to-gas ratio and mean 
extinction per unit gas, resulting in a steady growth in IRX. Higher mass galaxies had their 
periods of peak star formation in the past, and the exhaustion of their star forming gas 
supply (by whatever mechanism) leads to an IRX which falls with time. We reiterate that 
our determination of IRX is obtained directly from the Far UV-to-Far IR ratio (the latter 
from 12-24 fim luminosity and a bolometric correction), and is independent of the extinction 
law. 



We can model this using the classical exponential SFR models introduced by iTinsley 



( 119681 ) and a very simple extinction model. The SFR obeys: 

SFR ~ e-(*-*")/" = e^^ (8) 

where ^ is a scale-free "age" parameter. 

For simplicity, we characterize the extinction as if it occurs in a simple foreground 
absorbing slab of dust, and that its strength tracks the amount of gas responsible for star 
formation: 

^FUV ~ ZTjgas (9) 

where Z is the metallicity (we assume that gas-to-dust ratio scales accordingly), and T,gas 
is the gas surface density. Ignoring inclination induced anisotropics, we hav e IRX 



lo^Qio AApuv) _ 1), Let US further assume a Schmidt-Kennicutt scaling law (IKennicuttl 
1989h SFR ~ Ej„,. Here we note /? ~ 1.5. Then 



^gas ~ e-«/^ (10) 

In a closed-box enrichment model, metallicity grows as 

Z = ylrifi-^, = y ^ (11) 
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for an exponentially d ecaying SFR. The gas fraction is figas and y is the average yield 
( ISearle fc Sargentlll972l ). Thus we expect 



Afuv = Co^e-^/" 
Here Cq is a scaling constant appropriate for an average inclination. 



(12) 



This can be generalized to the leaky-box case (jHartwicklll976l ) in which the outflow is 
proportional to the star formation rate J^wind = —c SFR. In this case 



1 + c 



1 + c 



(13) 



If c < 1 then an accreting-box case (IBinney &: Merrifleldlll998l ) with infall proportional to 
SFR would obtain. We do not consider other accreting-box scenarios. 

Finally, we need to relate the age parameter to galaxy mass. We make a simple 
ansatz that the age scales as mass to a constant power a, that a single formation time t/ is 
appropriate, and that the past age is reduced by the relative elapsed time from formation: 



M 

Mo 



tiz) 



m-tf 



(14) 



The speciflc star formation rate (SSFR) is given for this model very simply by: 

SSFR=-- ^ 



re? - 1 



(15) 



where r is a function of mass for a coeval population: r = {A4/A4o) ° (13.47 — tj). 



We jointly flt the (IRX) and (SSFR) vs. mass and redshift with flve free parameters: 
Co, tf, Mq, a, and /?. We restrict the flts to the mass range 9.5 < \ogA4 < 11.5, over 
which the survey appears reasonably complete at all redshifts (cf. below). There are 46 
independent data points and four free parameters. The result is a reasonable flt, with 
= 31, with Co = 4.5, tj = 2.3 Gyr, logA^o = 10.58, a = 0.72, and /? = 2.1. Note that 
for Schmidt-Kennicutt star formation law we expected (3 ~ 1.5. The flts are displayed in 
Figure [12] and [121 Using this we predict that the mass metallicity relation shifts towards 
higher masses r oughl y A(logA^) ~ 0.53 at z = 0.7 and 0.75 at z = 1, not inconsistent with 
Savaglio et al.l (120051 ). This also predicts a peak equivalent extinction at ^peak = P = 2.2 



^We label galaxies by their stellar mass. Over the redshift range we consider, a constant SFR would 
increase the stellar mass by 0.5 dex, or one mass bin. We ignore this subtlety in order to make the arithmetic 
simple for this very basic model. 



-17- 



of AFuv,peak — 3.0, OT {I RX) peak=^ ■'^ ■ The oiily Way to distinguish between the closed and 
open box cases is to provide an independent measurement or prediction of the proportionahty 
constant Cq, which is beyond the scope of this paper. 

We show the tight relationship between the age parameter ^, derived from the stellar 
mass and the best fit parameters: 

M ft{z) -2.37 



and the mean IRX at in Figure [TH There is an equally tight relationship between the 
"b-parameter" which in terms of our simple parameterization is: 

b =< SSFR > [t{z) - t(0)] = (17) 

This is displayed in Figure [151 Finally, we show the co-evolution of (SSFR) and (IRX) in 
Figure [161 

We can also define a "turnoff" mass where ^{Ait,z) = 1.0 so that 

Mtiz) = 10^°-^ ( ^['l ~ ^-^l ) (18) 
^ ^ \t{0) - 2.37 J ^ ' 

which rises from log A^j = 10.5 at z=0 to log = 11.24 at z=l. We note that \ogMt = 10.5 
is exactly the transition mass between star forming and passively evolving galaxies noted by 



Kauffmann et al.l (120031 ) in the SDSS spectroscopic sample. 



Metallicity and age are simply related in this picture (cf. equation [TTj) . [Johnson et al. 



(12007) have shown that IR X correlates well with metallicity. The mass- metallicity relation- 



ship (iTremonti et al.[ [20041 ) is a result of the lower net astration in galaxies with a younger 
effective age ^. This simple picture does not predict the observed saturation at high mass. 
This could be a result of selection effects, since metallicity can only be measured using emis- 
sion lines in galaxies with high SSFR. It is plausible that transition galaxies with low emission 
line equivalent widths display a higher metallicity than the actively star forming sample at 
high mass. It could also indicate a second quenching mechanism in addition to simple gas 
exhaustion, or a more complex enrichment picture than simple closed box evolution. 



5.2. Star Formation History and Evolution of the Blue and Red Sequence 



We begin by determining the mass function by summing over SSFR in the bivariate 
mass-SSFR function. The distribution in each redshift bin is displayed in Figure [T71 The 
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error bars include Poisson and cosmic variance, the latter from ISomerville et al.l (l2004j ). 
Note that these combined errors will be highly correlated between mass bins at a given 
redshift. This correlation is not well represented by the plotted error bars, since there is 
large covariance between mass bins in a given redshift bin. Note also that there is evidence 
for incompleteness in the lowest mass bin (log = 9.75) for z=1.0. Mass bins log < 9.25 
were not used in the modelling of the previous section, because of obvious incompleteness. 



We model the mass function as an evolving Schechter function: 

M V ( M 



0(A<,2) =lnlO 4>,{z) 



exp 



(19) 



where A^*(z) = 7W*(0)(1 + z)'^ and 0*(-2) = 0^,(O)(1 + z^. We use a version of chi-square 
minimization that accounts for the cosmic variance and covariance between mass bins in a 



single redshift bin (e.g., iNewman fc David (120021 )) The best fit parameters and one sigma 
errors (obtained by marginalizing over other parameters) are 



0,(0) = 0.0040(-0.0023,+0.0036)Mpc~^logio^"^ 

logA^*(0) = 10.86(-0.43,+0.33) 

a = -0.93(-0.10,+0.11) 

P = 2.4(-1.0,+1.2) 

7 = -1.5(-l.l,+0.5) 



(20) 
(21) 
(22) 
(23) 
(24) 



The non-zero values of (3 and 7 indicate evolution in the characteristic mass and the number 
density, although cosmic variance prevents tight constraints on the evolutionary indices. 
Note we have not parameterized an evolving low mass slope. 



W e can compare the results of this mass function analysis with the results of iBell et al. 
(120031 ) using a much larger cosmic volume at z ~ 0.15. Converted to our cosmology and 



corrected upward to a non-diet Salpeter IMF (simply dividing their result by 0.7), they 
obtain 0, = 0.0035(0.0020) Mpc'^ log^o and logA^, = 11.17(11.06) for all (late-type) 

galaxies vs. our values at z=0.125 0* = 0.0040(-0.0023, +0.0036) Mpc'^ log^f^M'^, and 
log A^,, = 10.86. The mass cutoff is in fair agreement, but our density is a factor of 2 larger 
if we compare to just their late-type (morphologically selected) sample. The mass cutoff, 
slope and density parameters are highly correlated, and indeed if we fix our mass cutoff 
(z= 0) at 11.1 ( which is within our errors) we find 0* = 0.0025 Mpc~^ log^g-^"^) close to 
the bell et al.l J2003h value. If we fix the slope to a = -1.1 we find 0* = 0.00 20 Mpc" ^ 
l ogin -M "^ and logAl* = 11.1. Thus our results at z ~ are consistent with iBell et al. 



(120031 ) within the errors quoted above (dominated by cosmic variance). 



We can also compare to the evolving iBorch et al.l (120061 ) mass function which uses all 
three COMBO-17 fields with a total of 25,000 galaxies. We correct their masses to our 
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standard Salpeter IMF by multiplying by 1.8, as they suggest. At z=0.9, they obtain </)*(z = 
0.9) = 0.0012(0.0005) Mpc'^ logioA^"^ logM^ = 11.08(11.00) and , again for all (blue 
color-selected) galaxies. Using our evolutionary parameters, we find (f)^{z = 0. 9) = 0.0015 



Mpc ^ log^gA^ ^ and logJ^^{z = 0.9) = 11.54. Our results are compared to iBorch et al. 



(120061 ) in Figure [T71 Surprisingly, our results agree well with theirs for the entire galaxy 
sample in all but the highest redshift bin. Our bin extends over (0.8jzjl.2), while theirs 
is (0.8izil.O). Our mass function shows a distinctly higher mass cutoff in the two highest 
redshift bins, which leads to the stronger evolution in the characteristic mass. 

Searching for an explanation of this difference, we note that a significant fraction of 
our sample has red colors. In the lowest redshift bin this is true even after the extinction 
correction. For example, if we eliminate galaxies with {SSFR)< 10^^'^'^, corresponding to 
corrected {NUV — H)q ab > 4.25, which from Figure [3] can be seen to exclude the tail of 
the distribution, the mass function density parameter at z=0 falls by a factor of two. This 
suggests that a significant fraction of the mass (and even the star formation rate) at low 
redshift is locked in high mass galaxies with low specific star formation rates and re d intrinsic 



(unextincted) colors. Some of these galaxies could be classified as early-type in iBell et al. 



(120031 ) because their observed colors are even redder (as they are massive and exhibit high 



IRX) while their morphologies could be dominated by an evolved, bulge-like c omponent. 



The steepness of the slope parameter obtained for color-selected blue samples ( IBell et al. 



20031 ) could also be a refiection of excluding extincted, reddened higher mass galaxies and 



transition galaxies that still show some star formation. 

The characteristic mass increases with redshift, with AlogTW = 0.7 ± 0.4, which can 
be compared to the model change of logA^t from 10.5 to 11.24. Thus the blue sequence 
mass function shows directly the effects of downsizing. The total mass of the blue sequence 
obtained by integrating the mass function is a declining function of redsh ift (since > —7) 



Pm ^ 4 . 8 X 1 0^ 2.8 X lO^M^ Mpc'^ (from z=0.8 to z=0). Others (IBorch et all I2OO6 



BlantonI l2006l ) have found that the blue-sequence mass is constant with time, and giving 
the uncertainty our results are not inconsistent with this conclusion. In our model this 
occurs because the mass of the blue sequence is moving from fewer massive galaxie s to larger 



numb e rs of lower mass ga laxies. The decrease in blue luminosity density (e.g., iBell et al. 



(I2OO4J ) : iFaber et al.l (120051 )) is a result of an increase in mass-to-light ratio due in turn to the 
falloff of the specific star formation rate. The total stellar mass appears to remain constant 
or decline in spite of the formation of new stars. Before exploring this point, we estimate 
the star formation rate density and its evolution. 



In order to calculate the star formation rate history while minimizing the effects of 
cosmic variance, we renormalize the observed mass function to the fit given in Figure [17] 
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and Equation [20] for all redshift bins. The results in each mass bin and the total are shown 
in Figure [181 We also calculate the star formation rate density evolution for the simple 
exponential model of the previous section. This is shown in the right-hand panel of Figure 
[THl Apparently the dominant galaxy population responsible for the fall in SFR density from 
z=l to z=0 is 10.5 < logA^ < 11.5 (cyan and green points). This is consistent with the 
conclusion that the characteris tic mass Mc derived in the previous section evolves from lO^'''^ 
to 10^^'^^. IZheng et al.l (120071). also using COMBO-17 and Spitzer data, and the stacking 
technique of I Zheng et al.l (l2006l ). have reached conclusions that in many ways are similar 



to ours, but differ in some important details. In particular, they find that the SFR in the 
highest mass bin (log M<11.25 when converted to our IMF) does not fall more steeply than 
that in the lower mass bins. 

Stars formed over < 2; < 1 will increase the total mass of the blue sequence unless 
these galaxies transition to the red sequence. In order to maintain the constant or declining 
blue-sequence mass, we calculate that the average mass flux over < 2; < 1 must be Pb^r — 
0.03 — 0.05 M.Q yr~^ Mpc~^. This is in agreement with the value we obtained by examining 
transition galaxies at ^ ~ 0.1 ([Martin et al.ll2007al ). 



6. Summary 

We have used COMBO-17, Spitzer, and GALEX data to study the co-evolution of the 
IRX and star formation for galaxies over the mass range of 9.5 < logM < 12.0 and the 
redshift range < 2; < 1.2. We have reached a number of interesting conclusions: 

1. IRX grows with stellar mass in a way that mirrors the mass-metallicity relationship. 
The rise of IRX with mass saturates at a characteristic mass, above which it appears 
to fall. 

2. The SSFR is is roughly constant up to the same characteristic mass, above which it 
falls steeply. 

3. The characteristic mass grows with redshift. 

4. At a given mass below the characteristic mass, the IRX grows with redshift. 

5. The mass and evolutionary trends of the IRX and SSFR are reasonably fit by a simple 
gas-exhaustion model in which IRX is determined by gas surface density and metallic- 
ity, metallicity grows with time following a closed box model, and SFR is determined 
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by the exponentially falling gas density. The SFR time constant scales with the mass 
as r ~ M-°-^. 

6. The characteristic mass is a "turnoff" mass indicating galaxies that are starting to 
move off the blue sequence. 

7. The mass-metallicity relationship is understood to be determined largely by he char- 
acteristic age of the galaxies. The mass-IRX relationship is also influenced by gas 
exhaustion above the turnoff mass. 

8. The factor of 6-8 rise in SFR density to z=l is predominantly due to galaxies in the 
mass range 10.5 < logM < 11.5, the turnoff mass over the < z < 1 redshift range. 

These observations show directly the steady build-up of heavy elements in the interstellar 
media of evolving galaxies, and that the infrared excess IRX represents an excellent tool for 
selecting similar mass galaxies at various stages of evolution. In particular, galaxies at early 
stages in their evolution can be selected by their low IRX (Figure [Ml) . 

It is important to stress that these trends were uncovered by considering the average 
properties, notably the IRX, of large numbers of galaxies. A more sophisticated treatment 
would study the detailed distribution of properties, for example the spread in SSFR in a 
given mass bin. The simple scaling model of course predicts no spread at all for a coeval 
population. This distribution may yield information about the burst and formation history 
of galaxies. For example it will be very interesting to study galaxies with an unusually low 
IRX at a given epoch and mass to determine whether they are more recently formed. It 
will also be interesting to compare this very simple picture with the results of semianalytic 
models combined with the latest numerical simulations. Finally, it is critical to improve 
the observational basis of this work, most notably with a better understanding of the FIR 
bolometric correction and its evolution, with a larger and deeper sample of galaxies, and 
by extending the redshift range of this approach to determine whether this simple picture 
continues to apply during the major epoch of star formation. 

GALEX (Galaxy Evolution Explorer) is a NASA Small Explorer, launched in April 
2003. We gratefully acknowledge NASA's support for construction, operation, and science 
analysis for the GALEX mission, developed in cooperation with the Centre National d'Etudes 
Spatiales of France and the Korean Ministry of Science and Technology. 

FaciUUes: GALEX, SDSS 
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0.05<z<0.2 IRX 0.2<z<0.4 IRX 




0.4<z<0.6 IRX 0.6<z<0.8 IRX 





Fig. 1. — Volume corrected bivariate color-magnitude distribution 0(M/^, NUV — H). Con- 
tours are equally spaced in log0, with 10 divisions from —4 < log0 < —2. Colors give IRX, 
which is the log of the FUV-to-FIR luminosity ratio. Hashed color bins are upper limits to 
the IRX. Dotted lines give loci of constant M7V(7y=-15. 5,-18. 0,-20. 5. Dashed diagonal line 
indicates rriNuv = 26 limit in center of redshift bin. 
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Fig. 2. — Volume corrected bivariate extinction corrected color-magnitude distribution 
(p^Mn.o, {NUV — H)o). Contours are equally spaced in log0[Mpc~'^], with 10 divisions 
from —4 < \og(j) < —2. Colors give IRX, which is the log of the FUV-to-FIR luminos- 
ity ratio. Hashed color bins are upper limits to the IRX. Dotted lines give loci of constant 
Mjvc/v,o=-15.5,-18.0,-20.5. 
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NUV-H (NUV-H)(, 

Fig. 3. — Distribution of observed NUV-H (LEFT) and extinction-corrected (NUV-H)o 
(RIGHT). Color gives redshift bin: 0.05 < z < 0.2 (purple), 0.2 < z < 0.4 (blue), 
O.A< z < 0.6 (cyan), 0.6 < 2: < 0.8 (green), and 0.8 < z < 1.2 (red). 
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Fig. 4. — LEFT: Mass-to-liglit ratio at H-band vs ( NUV-H) for expon e ntially declining 
star formation histories e~'^*, obtained using model of iBruzual fc Charlotl (120031 ). RIGHT: 
Specific star formation rate (SSFR) for exponentially declining star formation histories. 
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Fig. 5. — Volume corrected bivariate Mass-SSFR distribution ^{M^SSFK). Contours are 
equally spaced in log0 [Mpc~'^], with 10 divisions from —4 < log0 < — 1. Colors give IRX, 
which is the log of the FUV-to-FIR luminosity ratio. 
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Fig. 6. — Volume corrected bivariate Mass-SSFR distribution weighted by the SFR, 
SFR*(f){M, SSFR). Contours are equally spaced in log^Fi? * [Mq yr'^ Mpc'^], with 10 
divisions from —4 < log SFR * (p < —1. Colors give IRX, which is the log of the FUV-to-FIR 
luminosity ratio. 
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Fig. 7. — Average IRX (IRX) (linear version plotted on a logarithmic scale) vs. mass in 
each redshift bin. LEFT: numerical average. RIGHT: SFR-weighted average. Errors are 
derived from bootstrapping. 
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Fig. 8. — Average IRX (IRX) (linear version plotted on a logarithmic scale) vs. redshift 
in each mass bin. LEFT: numerical average. RIGHT: SFR-weighted average. Errors are 
derived from bootstrapping. 
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Fig. 9. — Average specific star formation rate {{SSFR)) in eacli mass and redsliift bin. 
Errors are derived from bootstrapping. 
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Fig. 10. — Average IRX {{IRX)) vs. redshift in each mass bin using various techniques 
as indicated by hue styles. SOLID: 1. Basehne: NUV<26.0, r<24.0 DOTTED: Case 
2: NUV<25.0, r<24.0. DASHED: Case 3: NUV<27.0, r<24.0. DASH-DOT: Case 4: 
NUV<26.0, r<25.0 DASH-DOT-DOT: Case 5: all r<24.0 objects, whether or not detected 
in NUV. 
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Fig. 11.— 



Same as Figure 5, with color-coding changed to axis ratio b/a. 
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Fig. 12. — Simple SSFR-IRX model fits (lines) to (IRX) (linear version plotted on a 
logarithmic scale) vs. z (LEFT) and (IRX) vs. mass (RIGHT). LEFT: Color gives mass, 
with 9.5 < logM < 10.0 (purple), 10.0 < logM < 10.5 (blue), 10.5 < logM < 11.0 (cyan), 
11.0 < logM < 11.5 (green), 11.5 < logM < 12.0 (red). RIGHT: Color gives redshift bin: 
0.05 < 2 < 0.2 (purple), 0.2 < z < 0.4 (blue), 0.4 < ^ < 0.6 (cyan), 0.6 < ;z < 0.8 (green), 
and 0.8 < z < 1.2 (red). 
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Fig. 13.— Simple SSFR-IRX model fits (lines) to (SSFR) vs. z (LEFT) and {SSFR) vs. 
mass (RIGHT). LEFT: Color gives mass, with 9.5 < logM < 10.0 (purple), 10.0 < logM < 
10.5 (blue), 10.5 < logM < 11.0 (cyan), 11.0 < logAf < 11.5 (green), 11.5 < logM < 12.0 
(red). RIGHT: Color gives redshift bin: 0.05 < z < 0.2 (purple), 0.2 < z < 0.4 (blue), 
O.A< z < 0.6 (cyan), 0.6 < 2; < 0.8 (green), and 0.8 < z < 1.2 (red). 
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Fig. 14. — IRX (linear version plotted on a logarithmic scale) vs. galaxy age parameter ^ for 
all masses 9.5 < logA^ < 12.0 and all redsliifts. Color gives mass, with 9.5 < logM < 10.0 
(purple), 10.0 < logM < 10.5 (blue), 10.5 < logM < 11.0 (cyan), 11.0 < logM < 11.5 
(green), 11.5 < logM < 12.0 (red). Symbol type gives redshift: 0.05 < z < 0.2 (plus), 
0.2 < 2 < 0.4 (star), 0.4 < ;z < 0.6 (diamond), O.Q < z < 0.8 (triangle), and 0.8 < ^ < 1.2 
(square) . Solid line is from equation |4] and [121 
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Fig. 15. — Average b-parameter {{SSFR) times time since formation to) vs. galaxy age 
parameter ^ for all masses 9.5 < logAI < 12.0 and all redshifts. Color gives mass, with 
9.5 < logM < 10.0 (purple), 10.0 < logM < 10.5 (blue), 10.5 < logM < 11.0 (cyan), 
11.0 < logM < 11.5 (green), 11.5 < logM < 12.0 (red). Symbol type gives redsliift: 
0.05 < ;z < 0.2 (plus), 0.2 < ^ < 0.4 (star), OA < z < 0.6 (diamond), 0.6 < z < 0.8 
(triangle), and 0.8 < z < 1.2 (square). Solid line is from equation [T5l 
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Fig. 16. — IRX (linear version plotted on a logarithmic scale) vs. b-parameter 
{{SSFR) times time since formation to) for all masses 9.5 < logTW < 12.0 and all red- 
shifts. Color gives mass, with 9.5 < logM < 10.0 (purple), 10.0 < logM < 10.5 (blue), 
10.5 < logM < 11.0 (cyan), 11.0 < logM < 11.5 (green), 11.5 < logM < 12.0 (red). Sym- 
bol type gives redshift: 0.05 < z < 0.2 (plus), 0.2 < z < 0.4 (star), 0.4 < z < 0.6 (diamond), 
0.6 < z < 0.8 (triangle), and 0.8 < z < 1.2 (square). Solid line is from equation H], [12], and 
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Fig. 17. — Mass distribution function (/)(logA^) vs. logA^ in each redshift bin. Points and 
error bars give mass function derived in this paper. Sohd hne gives res ult of fit to evolvin g 
Schechter function. Colored lines give Schechter functions derived by iBorch et al.l (120061 ). 
Green (dotted) is for entire galaxy sample, red for red-sequen c e gala xies, and blue for blue- 
cloud galaxies. Note that the z=0 point is from iBell et al.l (120031 ). and our z~1.0 bin is 
compared to the iBorch et al.l (120061 ) z=0.9 bin. 
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Fig. 18. — Star formation rate density [A4q yr~^ Mpc~^] vs. redshift derived assuming a 
non-evolving mass function. Error bars are based on bootstrap but do not include cosmic 
variance, which adds a relative error of 0.37 to all the measurements. LEFT: density in each 
mass bin, with 9.5 < logM < 10.0 (purple), 10.0 < logM < 10.5 (blue), 10.5 < logM < 11.0 
(cyan), 11.0 < logM < 11.5 (green), 11.5 < logM < 12.0 (red), and total (black). RIGHT: 
Compared with simple IRX-SSFR model, plotted in each mass bin (colors as in left panel) 
and total (black). 
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Table 1. Source Detection and Matching Statistics 



Bands Number 



R<24 15882 

IRACl > 0.0005 mjy 13754 

MIPS24 > 0.02 mJy 3098 

R<24, NUV<26.0 10298 

R<24, NUV<26.0, FUV<26.0 4356 

R<24, IRAC1>0.0005 8196 

R<24, MIPS24>0.02 2090 

R<24, NUV<26.0, MIPS24>0.02 1481 



R<24, NUV<26.0, IRAC1>0.0005, MIPS24>0.02 1274 
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Table 2. log W ^a^ + bx log{L[X]/10^°LQ) 



X 


ax 


bA 


12 


0.710 


0.037 


13 


0.705 


0.036 


14 


0.778 


0.033 


15 


0.860 


0.030 


16 


0.884 


0.026 


17 


0.878 


0.023 


18 


0.874 


0.020 


19 


0.880 


0.016 


20 


0.874 


0.011 


21 


0.873 


0.007 


22 


0.863 


0.003 


23 


0.854 


0.000 


24 


0.824 


-0.004 



